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Abstract 



We present results of a search for R-parity-violating decay of the neutrahno 
Xi, taken to be the Lightest Supersymmetric Particle. It is assumed that this 
decay proceeds through one of the lepton-number violating couplings Aj^^ 
(j = 1,2; k = 1,2,3), and that R-parity is conserved in all other production 
and decay processes in the event. This scenario provides two muons and 
four jets in the final state. This search is based on 77.5 it 3.9 pb~^ of data, 
collected by the D0 experiment at the Fermilab Tevatron in pp collisions 
at y^=1.8 TeV in 1992-1995. Background expected from standard model 
processes amounts to 0.18 ± 0.03 ± 0.02 events. In the absence of candidate 
events, the result is interpreted in terms of limits on squark and gluino masses 
within the framework of the minimal low-energy supergravity supersymmetry 
model. 

A search for events with multiple leptons and jets is an effective way to look for new physics 
because such events do not suffer from large standard model (SM) backgrounds. These 
events can provide evidence of R-parity-vioIating (RPV) decays of supersymmetric (SUSY) 
particles [|l],||. R-parity is a discrete multiplicative quantum number that distinguishes 
SM particles from their SUSY partners. It is defined as R = (_ i)3-B+l+2S^ where B, L, 
and S are the baryon, lepton, and spin quantum numbers, respectively. R is +1 for SM 
particles and —1 for the corresponding SUSY particles. Originally, conservation of R-parity 
was imposed on supersymmetric theories because the combination of lepton-number and 
baryon-number violating couplings in the Lagrangian could have generated several rare or 
forbidden processes at unacceptably high rates. One such example is the decay of the proton. 
However, rapid proton decay as well as other rare decays can be prevented by not allowing 
simultaneous violations of baryon and lepton numbers. Thus, a small violation of R-parity 
cannot be excluded. 

The Yukawa coupling terms in the superpotential that induce R-parity violation are: 

XijkLiLjEk + \ji^LiQjDk + XijiJJiDjDk, 

where L and Q are the SU(2)-doublet lepton and quark superfields; E, U, and D are the 
singlet lepton, up-type quark, and down- type quark superfields, respectively; and i, j, and 
k are the generation indices. Since A and A are antisymmetric in the first two and last 
two indices, respectively, there are in total 45 possible couplings. For experimental searches 
it is usually assumed that only one of the 45 couplings is non-zero. Since experimental 
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upper bounds on these couplings from low-energy measurements are quite stringent 0], it 
is further assumed that R-parity violation manifests itself only in the decay of the lightest 
super symmetric particle (LSP). At the same time, these couplings are assumed to be strong 
enough so that the LSP is unstable and decays within the detector, close to the interaction 
vertex, which sets the scale for A at ^ 10^'^. A previous study at D0 |^ in the dielectron 
+ jets channel, searched for such a decay for non-vanishing A^^ (j = 1,2 and k = 1,2,3) 
couplings in the framework of the minimal low-energy supergravity supersymmetry model 
(mSUGRA) 0, with as the LSP. This model contains five parameters: a common mass 
for scalars (mo), a common mass for gauginos {mi/2), a common trilinear coupling [Aq, 
specified at the grand unification scale), the ratio of the vacuum expectation values of the 
two Higgs doublets (tan/3), and the sign of the Higgsino mass parameter (/i). The LSP 
decay to a charged lepton and two quark jets involving one of the X^j^. couplings is a viable 
mode for searching for SUSY at the Tevatron for the following reasons. The LSP can be 
produced either directly or through cascade decays from squarks or gluinos and can subse- 
quently decay into a lepton and two quarks. The branching fraction of this decay depends on 
the composition of the LSP, which in turn depends on the mSUGRA parameters described 
earlier. Studies have shown that at the energy of the Tevatron, the amount of signal in any 
of the lepton + jets decay channels of the LSP can be substantial for a large range of values 
of the mSUGRA parameters Also, such events will not contain any missing energy, 

thus making it easier to search for a RPV signal. We report a study similar to the previous 
one m, for finite A222 coupling (the study is equally valid for all the Aaj^ couplings with 
j = 1,2 and k = 1,2,3), based on a signature of two energetic muons and four energetic 
jets. There are several standard model processes that mimic this signature, e.g., 7*/2' — /z/i, 
Z ^ TT ^ /i/i, ti —>■ fifi, WW — >• /i/x, and accompanying jets. 

The D0 detector has been described elsewhere p|. The most important parts for this 
analysis are the uranium/liquid- argon calorimeter and the muon system. A cone algorithm 
with a cone radius of 0.5 in the rj-(j) space, where r] is the pseudorapidity and (j) is the az- 
imuthal angle, is used for jet identification 0. Muons are identified as tracks that leave 
minimum ionizing energy in the calorimeter, and are reconstructed in the muon system. An 
integrated luminosity of 77.5 ± 3.9 pb^^ collected with the D0 detector during the 1992- 
1995 Tevatron run at ^/s = 1.8 TeV is used for this analysis. The data are required to 
satisfy a trigger demanding one muon {px > 10 GeV/c, \ri\ < 1.7), and one jet {Ex > 15 
GeV, |?7| < 2.5). In the offline analysis, an event is selected only if it has at least two muons 
within \ri\ < 1.7 {pt > 15 GeV/c for the first muon, and pr > 10 GeV/c for the second 
muon), and at least four jets within \ri\ < 2.5 and with Et > 15 GeV. The muons and jets 
are required to satisfy standard D0 selection criteria [|TU|,|n|. The muons are also required 
to be isolated from jets by a distance > 0.5 in the r/-^ plane (this rejects muons coming from 
heavy-flavor decays, pions decaying in flight, and pion-induced punchthroughs) . In addition, 
several other criteria are imposed to minimize background. The aplanarity [|lT| of the jets 
in each event is required to be greater than 0.03. The invariant mass of the two muons 
is required to be greater than 5 GeV/c^, which helps to reject low-energy resonances (e.g. 
J/ip) and spurious combinations of muon tracks. Hjn, the scalar sum of Et of all muons and 
jets that pass kinematic and flducial requirements, is required to be greater than 150 GeV. 

Of the original 230,688 events passing the trigger requirements, none survive the above 
selections. The expected backgrounds from the two main SM channels, fifj,) + jets 
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and > /i/i) + jets, are shown in Table |, along with their statistical (first) and systematic 
(second) uncertainties. The contribution to background events from Z production is esti- 
mated from a sample of 21,000 Z + jets events, generated using VECBOS |T^. A total of 
254,000 events, generated with herwig 
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are used to estimate the contribution from 

which 



this background. The D0 detector is simulated using a GEANT-based package 
provides efficiencies of the selection criteria for signal and background events. We illustrate 
in Fig. |I|, the effect of one of the selection criteria (number of jets in an event) on events 
at a typical signal point (mo=140 GeV/c^, mi/2=90 GeV/c^, Ao=0, tan/?=2, /i < 0) and 
on events from the background channel Z{-^ fifi) + jets. The arrow in Fig. |I| indicates the 
minimum number of jets in accepted events. The instrumental background, which arises 

TABLE I. Summary of major backgrounds. First error is statistical and second error is sys- 
tematic. 



Background 


Expected events for 77.5 pb ^ 


process 




fifj,) + jets 


0.140 ±0.031 ±0.015 


tt{-^ fifi) + jets 


0.042 ±0.002 ±0.013 


Total 


0.182 ±0.031 ±0.020 



from misidentification of jets as muons, is negligible in this analysis. As can be seen from 
Table |I|, the expected number of background events is quite small. The statistical error arises 
from a combination of fluctuations in the Monte Carlo events and uncertainties in the muon 
and jet identification efficiencies. The systematic error arises due to uncertainty in the jet 
energy scale and in the values of production cross sections. 

Signal events are generated with isajet , modified to incorporate RP V decays based 
on the formalism of Ref. 0. For each signal sample, the value of efficiency multiplied by 
the branching fraction of — > > two muons and > four jets is estimated in the same way 
as described above for the SM background. Table || shows these values and the event yields 
expected from an integrated luminosity of 77.5 pb~^ for several points in the (mo, mi/2) 



TABLE II. Efficiency (e) multiplied by branching fraction {B), and expected event yield {N), 
for several points in the (mo, parameter space (for tan/?=2, Aq=0, and /i < 0). 



mo (GeV/c^) (GeV/c^) eB{%) {N) 

100 0.60 ± 0.07t^:^^ 3.0 ±0.4 

80 90 0.74 ± 0.08t[]:[]^ 2.7 ±0.3 

80 110 0.34 ± 0.04l°;°3 0.6 ±0.1 

190 90 0.78 ± 0.06l[]:[|| 2.1 ±0.2 

260 70 0.42 ± 0.04l[];[j^ 2.7 ±0.3 

400 90 0.31 ± 0.04l[];[j^ 0.8 ±0.1 
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parameter space. 

Since the expected SM background is compatible with absence of observed events, we 



c 

<U 300 
> 



5 6 7 

No. of Jets 



FIG. 1. Distribution of the number of jets per event at a typical signal point (mo=140 GeV/c^, 
mi/2=90 GeV/c^) (dashed line) and background sample, > /x^u) + jets (solid line) for 77.5 pb~^ 
integrated luminosity. The vertical arrow indicates the position of the applied cut. 



proceed to determine the region in mSUGRA space that can be excluded. An upper limit 
at the 95% confidence level (CL.) on the cross section for signal is obtained for each point 
in the (mo, 1711/2) plane for fixed values of Aq=0, fj, < 0, and tan/? = 2 and 6. A technique 
based on Bayesian statistics is used for this purpose, with a fiat prior for the signal cross 
section and Gaussian priors for luminosity, efficiency, and expected background. The hmits 
on the measured cross section are then compared with the leading-order SUSY prediction 
given by isajet, to find an excluded region in the (mo, m.1/2) plane. Figs. and | show the 
regions of parameter space (below the bold lines) excluded at the 95% CL. for tan/3 = 2 
and 6, respectively. 

The shaded areas in the left-hand corners of the figures indicate the regions where either 
the model does not produce electroweak symmetry breaking or the lightest neutralino is not 
the LSP. The area in the (mo, mi/2) plane excluded by experimental searches at LEP 
already extends beyond the shaded areas. The exclusion contour in Fig. ^ follows essentially 
a contour of constant squark mass (m^ = 260 GeV/c^) for low mo values. This is because 



pair production of squarks is the dominant SUSY process that contributes to the signal in 
that region. Production of gluinos, X2, and Xi becomes dominant at larger values of mo, 
where the masses and production cross sections of these particles are approximately inde- 
pendent of mo. The exclusion contour therefore becomes approximately independent of mo 
for mo > 250 GeV/c^. 

The value of Aq does not affect the results significantly, since it changes only the third 
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m„ (GeV/c^) 



FIG. 2. Exclusion contour in the (mo, mi/2) plane for tan/3 = 2, fj, < 0, Aq = 0, and finite Agj^ 
(j = 1,2; k = 1,2,3) coupling. The region below the bold line is excluded at the 95% C.L. The 
cross hatched region is excluded for theoretical reasons (see text), rriq and rrig denote the squark 
and gluino masses, respectively. 

generation sparticle masses. Both for /i > 0, and for higher values of tan/5 (see Fig. 3 for 
the exclusion contour at tan/5 = 6), the sensitivity of this search diminishes, because of the 
change in the composition of the LSP, which leads to a decrease of the branching fraction 
of the LSP into muons 0. 

In conclusion, we have searched for RPV decay of the neutralino Xi i^ito a muon and 
two jets in 77.5 pb~^ of data. No candidate events were found. This result is presented as 
an exclusion contour in the mSUGRA (mo, mi/2) parameter space for Ao=0, tan/3=2 and 
6, and /i < 0. In particular, for tan/5 = 2, squark masses below 240 GeV/c^ (for all gluino 
masses) and gluino masses below 224 GeV/c^ (for all squark masses) can be excluded. For 
equal masses of squarks and gluinos the mass limit is 265 GeV/c^. 

We thank the staffs at Fermilab and collaborating institutions, and acknowledge support 
from the Department of Energy and National Science Foundation (USA), Commissariat 
a L'Energie Atomique and CNRS/Institut National de Physique Nucleaire et de Physique 
des Particules (France), Ministry for Science and Technology and Ministry for Atomic En- 
ergy (Russia), CAPES and CNPq (Brazil), Departments of Atomic Energy and Science and 
Education (India), Colciencias (Colombia), CONACyT (Mexico), Ministry of Education and 
KOSEF (Korea), CONICET and UBACyT (Argentina), The Foundation for Fundamental 
Research on Matter (The Netherlands), PPARC (United Kingdom), Ministry of Education 
(Czech Republic), and the A. P. Sloan Foundation. 
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FIG. 3. Exclusion contour in the (mo, ^1/2) plane for tan/3 = 6, /U < 0, Aq = 0, and finite 
(j = 1, 2; A; = 1, 2, 3) coupling. 
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